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Abstract:  The  influence  of  a  S  adlayer  on  CO  adsorption  onto  Ni(100)  is  examined. 
Tight  binding  extended  Huckel  calculations  on  a  three  layer  model  slab  indicate  that 
the  interadsorbate  separation  distance  determines  not  only  the  mechanism,  but  also  the 
effect  of  the  interaction.  If  the  C-S  distance  is  short,  sulfur  induces  site  blockage  of  CO 
chemisorption  by  means  of  a  direct,  repulsive  interadsorbate  mechanism.  If  the  separation 
is  increased  beyond  the  normal  S-C  bond  range,  the  sulfur  adatoms  work  indirectly  via 
modification  of  the  electronic  structure  of  the  substrate.  This  is  a  form  of  through-bond 
coupling.  It  is  consistent  with  the  well-documented  sulfur  poisoning  of  CO  adsorption 
and  its  usual  explanation  via  relative  electronegativities  of  adsorbates,  but  there  are  some 


conceptual  differences.  At  longer  coadsorbate  separations,  there  is  an  interesting  reversal 
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The  effect  of  atomic  adsorbates  on  the  chemisorption  of  small  molecules  and  on  the 
rate  of  certain  catalytic  reactions  is  dramatic1.  Electronegative  impurities,  such  as  halo¬ 
gens  and  chalcogens,  tend  to  hinder  these  processes.  In  contrast,  alkali  metals  and  other 
electropositive  elements  can  function  as  promoters.  A  great  variety  of  experimental  and 
theoretical  studies  have  been  conducted  with  the  motivation  of  understanding  the  poison¬ 
ing  and  enhancement  effects  of  coadsorbates.  For  example,  the  adsorption  of  CO  onto 
sulfided  surfaces  is  often  used  as  a  model  for  the  sulfur  poisoning  of  Fisher-Tropsch  hydro¬ 
carbon  catalysis  from  CO  and  Ha.  We  will  focus  on  S/CO  coadsorbates  on  nickel  surfaces, 
as  the  body  of  experimental  work  dedicated  to  these  systems  is  voluminous. 

As  is  true  of  other  coadsorbate  systems,  the  basic  nature  of  the  interactions  between 
surface  species  is  a  matter  of  controversy.  Opinions  are  substantially  polarized  between 
two  extremes.  The  mechanism  is  described  either  as  dominated  by  delocalized  long-range 
effects  which  allow  a  single  impurity  adatom  to  modify  many  adsorption  sites1,  or  alter¬ 
natively  as  mediated  by  local  bonding  and  short-range  site  blockage9.  Using  the  army 
of  acronymical  methodologies  availabls  to  surface  science,  experimental  evidence  can  be 
found  to  support  either  side.  We  list  only  a  few  examples.  Goodman  and  Kiskinova2,4, 
and  Erley  and  Wagner*  advocate  the  long-range  theory  based  on  the  well-documented 
nonlinearity  of  both  the  CO  saturation  coverage  and  CO/Hj  catalytic  methanation  as  a 
function  of  preadsorbed  sulfur  on  Ni(100)  and  Ni(lll).  On  Ni(100),  a  one  fifth  sulfur 
monolayer  causes  an  order  of  magnitude  decreases  in  the  rate  of  methanation.  Goodman6 


concludes  that  each  sulfur  atom  affects  some  ten  Ni  surface  atoms.  Madix  et  al.  ,7*$  im¬ 
plicate  local  site  blockage  to  explain  similar  results.  Gland  et  al.B  favor  the  short-range 
model  for  Ni(100)  in  light  of  high  resolution  electon  energy  loss  spectroscopy  (HREELS) 
and  temperature  programed  desorption  (TPD)  studies.  The  infrared  reflection-adsorption 
spectroscopy  (IRAS)  work  of  Trenary  et  al.10  on  Ni(lll)  is  in  good  agreement  with  the 
latter. 

Theoretical  treatments  of  adatom  poisoning  and  promotion  of  CO  chemisorption  sore 
generally  presented  in  the  framework  of  the  Blyholder  model11.  The  adsorption  geometry 
is  widely  accepted  to  be  through  the  carbon  end,  exactly  or  nearly  perpendicular  to  either 
the  clean  or  preadsorbed  nickel  surface1**-*.  Regardless  of  the  specific  adsorption  site,  the 
chemisorptive  bond  in  the  Blyholder  model  results  from  electron  donation  from  the  CO 
5(7,  la,  into  the  empty  surface  levels  and  backdonation  from  the  surface  into  the  CO  2 x*, 
lb.  As  both  CO  levels  are  localized  on  the  carbon  end19,  they  are  at  least  geometrically 
well  suited  for  the  electron  transfers.  Sulfur  and  other  electronegative  adatoms  withdraw 


electron  density  from  the  surface.  Little  change  in  So  to  metal  donation,  2a,  is  presumed, 
but  metal  to  2k*  backbonding,  2b,  should  be  reduced,  and  hence  the  poisoning  effect.  As 
the  occupation  of  the  C-0  antibonding  2k*  decreases,  the  C-O  bond  strengthens  at  the 
expense  of  the  chemisorptive  bond.  Electropositive  adatoms  should  work  in  the  opposite 
manner.  Electron  density  is  pushed  into  the  surface,  thus  enhancing  its  backbonding 
capacity. 

Theoretical  studies  are  able  to  reproduce  poisoning  and  promotion,  and  most  indicate, 
at  minimum,  a  rudimentary  agreement  with  this  simple  MO  picture.  A  cornucopia  of 
calculation  methodologies  has  been  called  upon,  including  LC AO- type  on  clusters14*1'*  and 
monolayers18,  monolayer  SLAPW-type18,  cluster  LCGTO-Xa17"*,  the  effective-medium 
theory  applied  to  jellium  surfaces18-’*,  the  muffin-tin  approximation  applied  to  clusters19 
and  most  recently,  all-electron  local-density-functional  theory20.  Nevertheless,  there  is  no 
more  agreement  (as  to  the  range  or  detailed  nature  of  the  interactions)  between  the  various 
calculations  than  that  found  in  experimental  studies.  Benzinger  and  Madix18  suggest  that 
direct  interadsorbate  interactions  may  be  the  key  factor,  although  the  metal-mediated 
electron  transfers  predicted  from  the  Blyholder  mechanism  are  clearly  reproduced.  On  the 
other  hand,  Nsrskov  et  al.18-’*  find  long-range  electrostatic  forces  to  be  more  significant. 
The  FLAPW  results  of  Wimmer  et  al.20  straddle  the  contoversy.  A  long-range  interaction 
is  dominant  for  a  K  adatom,  but  for  S,  both  direct  and  substrate  mediated  effects  are 


surmised. 


In  this  contribution,  we  examine  the  S  +  CO  system  on  Ni(100)  with  the  extended 
Huckel  tight  binding  method1  la,i.  The  procedure  is  a  approximate  MO  method  with  well- 
known  deficiencies.  For  instance,  it  does  not  give  reliable  potential  energy  curves  nor  is 
it  able  to  properly  model  magnetism.  But  it  is  a  transparent  methodology  and  reveals 
clearly  the  basic  interactions  that  are  responsible  for  bonding.  The  calculations  we  report 
are  for  nonmagnetic,  spin  paired  metal  and  adsorbate. 

We  will  assemble  the  Ni/S/CO  system  from  its  various  component  parts.  Computa¬ 
tional  details  are  given  in  the  Appendix.  As  the  focus  of  previous  work  by  this  group  has 
been  the  adsorption  of  CO  onto  Ni(IOO)  and  other  surfaces”,  several  of  the  closely  related 
systems  —  the  bare  Ni(100)  surface,  the  unadsorbed  CO  net,  the  Ni(100)-c(2x2)CO  on-top 
and  Ni(lll)CO  bridging  —  are  well  characterized  and  our  discussion  will  be  appropriately 
brief. 

For  the  Ni(100)  substrate,  a  three  layer  slab  with  the  bulk  nickel  lattice  constant 
(a=3.524A”)  was  used.  The  nearest  neighbor  contact  is  2.49A  on  this  square  net  surface. 
The  calculated  density  of  states  (DOS)  exhibits  a  compact  d  band  between  —12  and 
— 8eV,  and  very  disperse  s  and  p  bands  between  —12  and  8eV”’14.  As  the  surface  atoms 
have  a  lower  coordination  number,  these  states  are  less  disperse  than  the  bulk  states. 
Consequently,  if  the  Fermi  level  falls  above  the  midpoint  of  the  d  block,  as  for  nickel,  the 
surface  is  negatively  charged  with  respect  to  the  bulk. 

Whether  the  three  layer  slab  is  an  appropriate  model  for  the  surface  can  be  estimated 
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from  the  charge  on  the  middle  layer.  An  infinitely  deep  system  will  have  a  bulk  charge 
approaching  zero.  With  a  55  k  point  set1',  we  compute  an  excess  charge  of  +0.180e~  per 
atom  for  three  layers,  +0.166e~  for  four  layers  and  +0.157e~  (middle  layer)  for  five.  As  the 
convergence  is  slow  in  this  range  of  reasonably  sized  unit  cells,  we  make  our  choice  on  the 
basis  of  computational  economics.  The  three  layer  substrate  was  used  for  all  calculations. 

Our  first  priority  in  the  study  of  the  S/CO  system  on  Ni(lOO)  is  to  establish  proper 
adsorption  geometries.  So  let  us  review  briefly  what  is  known  from  experiment.  The 
preference  for  atomic  adsorption  at  the  highest  coordination  site  is  nearly  universal  on 
clean  or  coabsorbed  transition  metal  surfaces”.  Specifically,  chalcogen  adsorption  at  the 
Ni(100)  four-fold  hollow  is  well  documented*7*"*.  CO  chemisorbs  molecularly  and  does 
not  dissociate  to  any  measurable  degree  on  late  transition  metal  surfaces2’.  On  clean 
Ni(100),  CO  site  preference  is  coverage  and  temperature  dependent29*"*.  At  saturation 
coverage  and  below  150K,  a  compression  structure  of  both  on-top  and  bridging  CO  is 
formed.  When  the  compression  structure  is  relaxed  to  a  0=\,  c(2x2)  monolayer  (or  if 
coverage  is  further  reduced),  only  the  on-top  mode  is  observed.  The  HREELS  vibrational 
spectrum9  is  characterized  by  a  single  C-O  stretching  frequency  at  2005cm"1.  If  sulfur  is 
added,  and  its  surface  concentration  is  increased,  this  dominant  peak  is  downshifted  and 
two  new  CO  modes  grow  in,  one  at  a  lower  frequency  and  the  other  at  a  higher  frequency. 
The  adsorption  geometry  of  the  new  modes  is  unknown.  Using  partial  thermal  desorption, 
the  high  frequency  peak  is  correlated  to  low  temperature  desorption,  the  low  frequency 


to  high  temperature  and  the  dominant  mode  to  an  intermediate  temperature.  The  CO 
layer  is  completely  desorbed  below  430K,  which  is  the  predominant  desorption  temperature 
from  the  clean  surface.  The  inverse  relationship  between  the  C-0  stretch  frequency  and 
desorption  temperature  is  consistent  with  the  greater  susceptibility  of  the  chemisorptive 
bond  strength  to  the  x,  rather  than  the  a  system,  as  advocated  in  the  Blyholder  mechanism. 
However,  we  expect  a  unidirectional  effect;  the  sulfur  should  reduce  x-backbonding.  The 
C-0  bond  strengthens  as  the  chemisorptive  bond  weakens,  thereby  creating  a  high  CO 
frequency,  low  desorption  temperature  mode.  Although  two  such  modes  are  observed,  a 
third,  apparently  originating  from  increased  x-backbonding,  does  as  well.  We  ambitiously 
hope  to  model  these  results  by  examining  the  effect  of  sulfur  coadsorption  on  the  binding 
of  CO  at  various  sites  and,  as  well,  gun  general  insight  into  the  mechanism  of  coadsorbate 
interactions. 

A  coverage  of  0= j,  p(2x2)  per  adsorbate  species  was  chosen.  It  is  in  this  coverage  range 
that  the  three  C-0  stretch  peaks  first  appear  simultaneously  in  the  HREELS  specta9.  Each 
species  is  laid  out  in  a  square  net,  lattice  constant  a=4.98A.  The  calculated  band  structures 
of  CO23  and  S  unsupported  nets  are  essentially  free  of  dispersion.  The  band  energies  are 
equivalent  to  the  molecular  orbital  energies,  thus  intraspecies  interactions  are  negligible  at 
this  distance.  To  help  us  see  the  important  orbitals  of  both  adsorbates  and  their  relative 
energies,  we  show  in  Figure  1  the  total  density  of  states  (DOS)  of  the  bare  Ni  slab,  and 
alongside,  the  important  frontier  orbitals  of  S  and  CO. 


One  of  the  prerequisites  for  the  electronegativity  arguments  concerning  poisoning  is 
that  the  electron  withdrawing  power  of  sulfur  be  observed  in  Ni(100)-p(2x2)S,  that  is, 
without  CO  coadsorption.  Ten  and  16  k  point  sets  were  used  for  substrate  plus  adsor¬ 
bate^)  systems  with  tetragonal  and  orthorhombic  symmetry  respectively  (see  Appendix 
for  calculations  testing  convergence).  The  sulfur  square  net  is  placed  on  the  substrate 
such  that  S  occupies  four-fold  hollow  sites  and  d(S-Ni)=2.19A30,  as  obtained  from  LEED. 
The  z  axis  is  consistently  oriented  perpendicular  to  the  Ni  plane,  and  the  x  and  y  axes 
parallel  to  the  surface  Ni-Ni  close  contacts.  The  charge  transfer  occurs  as  predicted;  the 
sulfur  adatoms  gain  0.756e~,  and  each  of  the  four  associated  surface  atoms  looses  0.355e~ 
relative  to  bare  Ni(100).  Although  the  charge  transfer  is  greatest  in  the  nickel  xz  and  yz 
levels  (-0.115e-  each)  and  the  sulfur  z  (-f0.481e-si),  the  projected  DOS  indicate  that  the 
strongest  interaction  is  with  the  sulfur  x  and  y  rather  than  z  levels.  The  p  band  of  the  S 
square  net  lies  at  —  13.2eV,  well  below  both  the  d  block  and  the  Fermi  energy  of  the  bare 
substrate  (— 8.7eV).  The  large  dispersion  provided  by  the  interactions  with  surface  xz  and 
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Figure  1.  In  the  left  panel,  the  bare  Ni(IOO)  DOS.  The  valence  orbitals  of  CO  and  S  are 
drawn  in  the  middle  and  right  panels  respectively. 
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yz  states,  such  as  the  representative  one  depicted  in  3,  forces  some  20%  of  the  S  x  and  y 
states  above  the  Fermi  energy  and  splits  the  major  S  x  and  y  peak  down  from  the  z  by 
~0.5eV.  Some  surface  xz  and  yz  character  appears  in  the  S  x  and  y  peak  near  —  14e V,  as 
seen  in  a  comparison  of  the  projected  DOS  integration  curves  in  Figure  2.  In  contrast,  the 
S  z  level  acts  essentially  as  an  inert  electron  sink. 


The  total  DOS  (solid  line  in  Figure  2)  appears  very  much  as  a  simple  overlay  of  the  S 
square  net  DOS  (represented  schematically  on  the  right  by  bars  as  the  median  energies) 
and  the  bare  Ni(100)  DOS*®.  The  binding  energies  of  the  sulfur  p  and  s  levels  agree  with 
those  obtained  by  UPS*7"  within  experimental  error.  The  S  adlayer  is  a  relatively  small 
pertubation  on  the  electronic  structure  of  the  Ni(100)  slab  as  a  whole,  but  does  cause 
substantial  charge  redistribution  and  rehybridization  at  a  local  level.  If  poisoning  does 
occur  via  the  substrate,  the  xz  and  yz  surface  levels  are  the  most  likely  mediators. 

We  consider  three  high  symmetry  sites  for  CO  coordination:  on-top,  bridging  and 
four-fold  hollow.  The  C-Ni  and  C-0  bond  lengths  are  determined  by  LEED  as  1.80 A  and 
1.15 A  respectively1*"  for  on-top  coordination,  the  only  one  known  to  exist  independently.  In 
general,  our  computational  method  cannot  be  trusted  to  predict  the  energetics  of  adsorbate 
site  preference.  So  we  must  make  some  geometric  choices.  For  the  bridging  and  4-fold  sites, 
all  calculations  on  both  clean  and  sulfur  coadsorbed  surfaces  were  performed  at  two  limiting 


Figure  3.  The  solid  line  indicates  the  Ni(100)-p(2x2)S  total  DOS.  The  dotted  and  dashed 
lines  are  the  integrals  of  the  surface  xx  +  yz  and  Sx  +  y  projected  DOS  respectively. 
The  S  unsupported  square  net  DOS  is  represented  by  the  median  energy  bars  on  the 


geometries:  C-Ni  nearest  neighbor  distance  of  1.80A,  and  C-Ni  surface  distance  of  1.80A. 
The  trends  are  identical  but  consistently  more  pronounced  for  the  first  option,  which  we 
choose  to  present  for  the  sake  of  clarity. 

Characteristic  of  the  rather  weak  chemisorption  of  CO  onto  nickel  and  other  late  tran¬ 
sition  metals,  the  CO  molecular  levels  for  the  observed  on-top  geometry  form  bands  which 
have  a  small  dispersion  and  are  at  nearly  the  same  energy  as  the  molecular  levels.  In  con¬ 
trast,  one  third  of  the  bridging  and  4-fold  2**  levels  are  pulled  into  the  d  block.  Specific 
features  of  the  DOS’s  will  be  discussed  later,  in  conjunction  with  the  coadsorbed  systems. 

Our  models  adhere  closely  to  the  Blyholder  model  of  cr-donation  and  x-backdonation. 
The  electron  density  shifts  from  bare  Ni(100)  or  molecular  CO  to  Ni(100)-p(2x2)CO  are 
presented  in  Table  1.  The  extent  of  5a  depopulation  falls  within  a  much  narrower  range 


than  does  the  amount  of  backdonation  into  2 The  a  interaction  contributes  to  the 
net  bonding  in  a  relatively  uniform  manner  for  all  sites.  On  the  other  hand,  the  C-0 
overlap  population  (o.p.)  follows  the  2x*  occupation  exactly  as  predicted.  Moving  through 
the  coordination  series  of  molecular,  on-top,  bridging  and  4-fold,  the  C-0  o.p.  falls  in 
conjunction  with  growing  2 w*  occupation.  The  CO  coordination  number  mediates  the 
electronic  shifts  at  the  associated  surface  atoms  the  largest  shift  coincides  with  the  lowest 
coordination.  The  average  electron  loss,  however,  correlates  well  with  the  gain  in  2rr*, 


Table  1.  Calculated  Results  for  Ni(100)-p(2x2)CO 


Atom(s)  or 
Orbitals 


CO 

1  Coordination  Mode 

on-top 

4-fold 

bridging 

Electron 

Density 

Changes^ 


CO  total 
CO  5<r 
CO  2**^ 

Ni,  coordinated 
N  i,  uncoordinated 
Ni,,  average 


+0.249 

+0.777 

-0.379 

-0.368 

+0.372 

+0.694 

-0.814 

-0.283 

-0.005 

— 

-0.200 

-0.283 

+0.638 

-0.398 


+0. 

-0.542 


-0.002 

-0.273 


Overlap 

c-o(*) 

1.044 

0.916 

Populations 

C-Ni, 

0.847 

0.359 

CO  binding  E  (eV)(d> 

+2.570 

+2.526 

+3.843 

•Relative  to  Ni(100)  or  molecular  CO.  Net  electron  gain  if  positive,  loss  if  negative. 


bOccupation  of  2tt*  and  other  degenerate  orbitals  are  given  for  each  individual  orbital. 
'Molecular  C-0  o.p. =1.208 

dE(Ni(100))  +  E (molecular  CO)  -  E(Ni(100)-p(2x2)CO) 
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but  is  independent  of  5a.  These  two  effects  underline  the  notion  that  the  chemisorptive 
interaction  is  governed  by  the  extent  of  x-backbonding  against  a  fairly  uniform  o-donation 
background.  Unfortunately,  we  cannot  compare  the  strengths  of  the  chemisorptive  bonds 
directly;  the  effect  of  coordination  number  on  overlap  population  is  generally  nonlinear. 

Atomic  electron  transfers  must,  however,  be  used  with  care.  In  this  case,  both  CO 
interactions  work  to  decrease  the  surface  electron  density.  Most  of  the  d  block  lies  just 
below  the  Fermi  energy  so  that  some  states  will  be  pushed  above  it  by  dispersion  from  any 
source,  whether  from  below  (5oj  or  above  (2x*)2*.  That  d  block  states  can  be  reused  above 
ff  by  interaction  with  2x*  merits  a  word  of  discussion.  Very  simplistically,  mixing  between 
a  fully  occupied  “d"  and  empty  “2jt‘"  states  will  create  a  “d+A2x*”  state  below  cf  and  a 
“2x*— Ad”  state  above  <r  (A  is  an  arbitrary  mixing  parameter).  The  net  charge  transfer 
will  be  out  of  “d”  and  into  Although  the  bare  surface  d  block  occupation  is  actually 

9.3e~  rather  than  10.0e~,  it  is  not  hard  to  imagine  that  simply  by  their  preponderance,  the 
occupied  states  will  dominate.  The  two  effects  are  most  cleanly  separated  at  the  on-top 
geometry.  The  orbitals  able  to  interact  with  the  5 a  do  not  interact  effectively  with  the  2 x*, 
and  vice  versa.  The  appropriate  surface  d  orbitals  (z*  for  5 a,  xs  and  yz  for  2x*)  experience 
nearly  identical  charge  shifts.  Since  the  5 a  occupation  is  very  similar  at  the  three  sites,  we 
can  safely  surmise  that  its  interaction  with  the  substate  is  constant  and  the  total  surface 
electron  loss  is  determined  by  the  2 x*. 

With  regard  to  the  electronegativity  arguments  concerning  sulfur  poisoning,  we  suggest 
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that  the  electron  withdrawing  power  of  CO  must  be  considered  as  well.  The  total  charge 


transfer  into  CO  is  approximately  one  third  of  that  into  the  S  of  Ni(100)-p(2x2)S  for  CO 


on-top,  but  nearly  equivalent  for  the  bridging  and  4-fold  coordinations.  Although  the  point 


must  be  reconsidered  for  the  coadsorbed  systems,  this  result  makes  rationalizations  based 


on  electronegativity  alone  less  self-evident. 


Four  Ni(100)-p(2x2)S+p(2x2)CO  systems  have  been  considered:  1)  CO  on-top,  4a,  2) 


CO  in  four-fold  hollow  retaining  tetragonal  symmetry  (4-fold  tetra),  4b,  3)  CO  in  four-fold 


hollow,  symmetry  reduced  to  orthorhombic  (4-fold  ortho),  4c,  and  4)  CO  bridging,  4d. 


0 

S — C 

i  y 


The  top  layer  of  the  three  layer  unit  cell  is  shown  in  the  schematics.  In  each  case,  we 


have  left  the  S  in  a  four-fold  hollow.  Behind  this  assumption  is  the  strong  propensity 


of  the  S  for  this  site;  to  our  knowledge,  no  deviation  from  such  adsorption  is  known 


experimentally.  CO,  on  the  other  hand,  is  known  to  have  low  barriers  to  moving  into 


alternative  coordination  model. 

The  evolution  of  the  C-O  and  C-Ni,  o.p.’s  (Ni1=surface  atom)  upon  sulfur  coadsorp¬ 
tion,  shown  in  Table  2,  are  intriguing  in  light  of  the  Blyholder  and  “anti-Blyholder”  modes 
observed  by  Gland9.  Just  as  predicted  by  the  Blyholder  mechanism,  9  coadsorption  in 

Table  2  here 

the  4-fold  ortho,  4c,  and  bridging,  4d,  systems  induces  C-O  bond  strengthening,  C- 
Ni,  bond  weakening  and  reduces  the  CO  binding  energy.  Although  the  absolute  binding 
energies  calculated  by  the  extended  Huckel  method  are  unreliable  and  we  are  unable  to 
predict  adsorption  site  preferences,  the  relative  binding  energies  at  the  same  site  can  be 
used  with  more  confidence.  We  tentatively  associate  the  4-fold  ortho  and  bridging  sites  to 
the  high  frequency  C-O  stretch,  low  temperature  desorption  mode.  Exactly  the  opposite 
effects  are  observed  for  the  4-fold  tetra,  4b.  This  site  corresponds  to  the  “anti-Blyholder” 
low  frequency  C-O  stretch,  high  desorption  temperature  mode.  That  the  interaction 
between  adsorbates  is  non-isotopic  is  not  without  precedence.  Other  workers33*’*  have 
proposed  that  the  interaction  is  an  oscillatory  function  of  the  separation  distance,  thus 
alternately  repulsive  or  attractive. 

Both  bonds  are  weakened  if  CO  is  bound  on-top,  and  the  CO  binding  energy  decrease 
so  sharply  that  the  Ni(100)-p(2x2)S+p(2x2)CO  system  is  less  stable  than  its  separated 
components.  The  anomalous  behavior  arises  directly  from  the  close  C-S  contact.  In 


16 


Table  2.  Overlap  Populations  k  Binding  Energies 

C-O  Overlap  Population  C-Ni,  Overlap  Population  ACO  binding 
CO  Site  Without  S  With  S  Without  S  With  S  Energy (eV)W 

on-top,  4a  1.044  0.929  0.847  0.579  +3.391 

4-fold  tetra,  4b  0.916  0.906  0.359  0.363  -0.050 

4-fold  ortho,  4c  0.916  0.928  0.359  0.336  +2.435 

bridging,  4d  0.942  0.974  0.631  0.581  +1.256 


this  geometry,  4a,  with  £  coverage  for  both  S  and  CO,  the  S-C(O)  separation  is  only 
1.83 A.  This  is  very  much  within  a  bonding  range;  indeed  typical  S-C  single  bonds  fall 
between  1.75  and  1.80A.  The  computed  C-S  o.p.  is  0.627,  which  is  the  largest  value  of 
any  site  (the  bridging  mode  is  a  distant  second  at  0.015).  The  surface  species  is  better 
described  as  molecular  SCO  than  as  coadsorbed  S  and  CO.  To  our  knowledge,  no  evidence 
of  SCO  formation  has  been  observed  for  any  coadsorbed  S/CO  surface.  The  energy  of 
this  coadsorbate  geometry  is  very  high.  We  believe  it  is  energetically  inaccesible,  and  that 
what  we  are  in  fact  modeling  is  local  site  blockage. 

Essentially  all  other  calculations  to  date  have  been  performed  at  similarly  short  C-S 
contacts.  This  distance  is  crucial;  at  1.83A  we  observe  site  blockage,  at  2.66A  (4-fold  ortho) 
and  2.78A  (bridging)  the  Blyholder  result,  and  at  3.64A  (4-fold  tetra)  an  “anti-Blyholder” 
effect.  The  smallest  distance  must  be  classified  as  short-range  and  the  large  S-C  o.p. 
stongly  indicates  a  direct  interadsorbate  effect;  the  largest  is  clearly  long-range  and  its 
corresponding  o.p.  is  negligible  (-0.003),  pointing  to  a  surface-mediated  interaction.  The 
question,  then,  may  not  be  what  is  the  absolute  range  of  interadsorbate  interactions  but 
rather  to  identify  the  type  of  interaction  active  in  each  range. 

Tibfo  3  htta 

An  important  point  in  our  analysis  is  that  the  2x*  occupation  follows  the  strength  of  the 
chemisorptive  bond.  The  trends  shown  in  Table  3  are  consistent  with  the  v-b&ckbonding 


Table  3.  Occupation  of  2* 


CO  Site 

Predicted  Changed 

Occupation 

Without  S 

Occupation 

with  s 

on-top 

— 

0.372 

0.515 

4-fold  tetra 

+ 

0.694 

0.717 

4-fold  ortho 

— 

0.694 

0.686 

bridging 

— 

0.598 

0.521 

‘Based  on  C-Ni,  overlap  population  trends. 


mechanism;  invariably  C-Ni,  o.p.  reduction  is  coincident  with  2r *  depopulation,  and  vice 
versa. 

The  relative  electron  withdrawing  power  of  the  two  adsorbates  is  reconsidered  in 
Table  4.  Relative  charge  transfer  is  defined  as  the  atomic  electron  density  in  Ni(100)- 
p(2x2)S+p(2x2)CO  minus  that  in  the  single  adsorbate  system  (Ni(lOO)-p(2x2)S  or  Ni(100)- 
p(2x2)CO).  A  positive  value  indicates  that  the  coadsorbed  species  aquires  electron  density, 
a  negative  value  indicates  a  net  loss.  According  to  this  indicator,  the  S  is  more  electroneg¬ 
ative  than  4-fold  ortho  or  bridging  CO,  but  more  electropositive  than  4-fold  tetra  CO.  The 
results  are  consistent  with  the  notion  that  a  more  electronegative  species,  the  4-fold  CO 
tetra,  will  accept  more  back  bonding. 

Table  i  hies 

Although  the  results  thus  far  are  consistent  for  each  site,  the  question  remains  as  to 
why  these  sites  behave  in  fundamentally  different  ways.  Let  us  begin  with  the  bridging 
geometry.  The  arguments  for  the  4-fold  ortho  site  are  much  the  same. 

The  composite  characteristics  of  the  Ni(100)-p(2x2)S+p(2x2)CO  bridging  total  DOS 
can  be  seen  clearly  in  Figure  3.  The  DOS  of  the  CO  and  S  unsupported  nets  are  represented 
by  bars  at  the  right.  In  spite  of  energy  shifts,  their  peakedness  is  substantially  retained;  a 
minimum  of  50%  of  the  adsorbate  levels  fall  in  their  major  peaks  as  labeled.  As  previously 
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Table  4.  Electron  Transfer  to  Coadsorbates  Relative  to  Single  Species  Systems^ 


*A  positive  value  indicates  an  electron  density  gain  relative  to  single  species  adsorption, 
and  a  negative  value,  an  electron  loss. 
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CO  2 7T* 


CO  5a 
3  3p 
CO  1 7T 

CO  4a 


S  3s 


Figure  3.  The  bridging  Ni(100)-p(2x2)S+p(2x2)CO  total  DOS;  major  peaks  are  labeled. 
The  ban  on  the  right  indicated  the  median  energies  of  the  CO  and  S  unsupported 
nets. 
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discussed,  the  direct  interaction  is  predictably  small  with  d(C-S)=2.78A  and  C-S  o.p. =0.015 
Some  mixing  occurs  between  the  CO  lx  and  the  sulfur  p  levels  since  they  are  energetically 
closely  matched.  But  the  lx  is  heavily  localized  at  the  oxygen  end  (see  Figure  1)  and  is 
therefore  relatively  inert  to  adsorption.  The  poisoning  must  occur  via  the  substrate. 

The  same  conclusions  can  be  drawn  from  the  electron  density  shifts,  of  which  a  selected 
few  are  listed  in  Table  5.  The  excess  charge  at  S  is  nearly  independent  of  CO  adsorption.  In 
contrast,  the  local  surface  charge  is  very  much  determined  by  the  adjacent  surface  species. 
With  only  one  exception,  every  atom  of  the  initially  negatively  charged  surface  (bare 


Ni(100))  loses  electron  density  as  adsorbates  are  added  sequentially.  The  bridged  atoms 
of  the  coadsorbed  system  suffer  most  with  respect  to  either  single  species  slab. 

Let  us  focus  on  the  xs  orbitals,  which  lie  along  the  bridging  direction.  As  discussed 
earlier,  the  electron  withdrawing  power  of  S  is  felt  most  strongly  at  this  orbital  and  its 
degenerate  partner  in  Ni(100)-p(2x2)S.  The  largest  perturbation  of  CO  adsorption  onto 
either  the  clean  or  sulfided  surface  occurs  here  as  well.  The  xz  is  geometrically  well  suited 
for  interaction  with  both  the  5<r  and  the  2x*x,  which  lies  along  the  bridging  direction, 
depending  on  the  phase  relationship  between  the  bridged  atoms.  The  out-of-phase  com¬ 
bination  can  participate  in  o-bonding,  5a,  and  the  in-phase  in  x-backbonding,  5b.  The 
interactions  are  energetically  favored  as  well.  The  out-of-phase  states  will  predominate  in 
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the  lower  half  of  the  d  band  (near  5 a)  since  they  are  bonding  between  the  metal  centers. 
Likewise,  the  in-phase  set  will  be  concentrated  in  the  upper  half,  adjacent  to  the  2**.  The 
interaction  with  the  2x*x  is  particularly  strong;  the  occupations  of  the  originally  degenerate 
molecular  orbitals  differ  by  0.297e~  and  the  median  energies  are  split  by  ~leV.  Later,  we 
will  use  the  phase  relationships  between  the  Ni  atoms  and  a  technique  we  call  COOP  to 
identify  the  a  and  ir  interacting  states. 

The  stepwise  evolution  of  the  x*  can  be  traced  in  Figure  4.  The  integral  of  the  projected 
DOS  is  compared  for  Ni(100),  Ni(100)-p(2x2)CO  and  Ni(100)-p(2x2)S+p(2x2)CO,  and 
overlayed  onto  the  total  DOS  of  the  latter.  The  initial  localization  of  the  states  entirely 
within  the  d  block  is  disrupted  by  both  o  and  ir  effects.  The  mixing  is  evident  from 
the  substantial  redistribution  of  levels  into  both  the  5 a  and  2 x*  areas.  Both  effects  are 
magnified  with  addition  of  the  S  adlayer. 

Figure  4  here 

If  the  integrals  of  the  Sc  projected  DOS  are  compared  for  the  surfaces  with  and  without 
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Energy  (eV) 


- Ni(100)-p(2x2)S+p(2x2)CO 

xz,  with  3  &  CO 
—  xz,  with  CO 
----  xz,  bare  surface 


CO  2n* 


CO  5 a 


Figure  4.  the  dashed,  dot-dashed  and  dotted  lines  are  respectively  the  integrals  of  the  xz 


projected  DOS  of  Ni(100),  Ni(100)-p(2x2)CO  and  Ni(100)-p(2x2)S+p(2x2)CO.  The 
solid  curve  is  the  total  DOS  of  the  latter.  The  median  energies  of  the  5 a  and  2? r* 
coadsorbed  levels  are  indicated  on  the  right. 


S,  we  can  see  that  sulfur  addition  pulls  some  of  the  5o  states  into  the  CO  1*  and  sulfur 
p  peaks,  but  does  not  substantially  alter  either  the  median  of  energy  or  the  distribution 
of  levels  near  «f.  On  the  other  hand,  a  similar  comparison  of  the  2n*x  integrals  reveals 
that  the  2x*x  level  is  strongly  distorted  and  the  median  is  pushed  up  ~leV.  We  must 
remember,  however,  that  these  effects  cannot  be  caused  by  a  direct  S-CO  interaction. 

How  does  the  depopulation  of  the  CO  2x*  come  about  on  the  sulfided  surface?  We 
saw  earlier  that  the  2 jt*  interacts  strongly  with  the  upper  part  of  the  metal  xz  band  (5b). 
The  S  p  orbitals  also  interact  with  the  same  part  of  the  band  and  push  it  up  (Figure  4). 
The  metal-2x*  interaction  becomes  stronger  since  the  energy  difference  is  reduced.  As  a 
result,  both  the  adsorbate  and  metal  levels  are  depopulated.  To  explain  why,  we  turn  to 
the  simple  two  level  interaction  diagram  used  so  commonly  for  molecular  systems,  and 
modify  it  to  indicate  the  inherent  level  dispersion  of  infinite  systems,  0.  The  schematic  is 


with  S 


upward  the  resultant  antibonding  levels,  which  contain  both  metal  and  2 it*  character. 
Contrary  to  molecular  systems  for  which  the  energy  of  the  HOMO  can  change  due  to 
interactions,  the  fermi  level  of  the  metal  surface  is  essentially  fixed,  regardless  of  the 
surface  species.  Our  three  layer  calculation  models  this  well;  the  greatest  shift  of  <f  after 
adsorption  is  0.13eV  from  the  value  of  the  bare  surface.  The  placement  of  Cf  in  6  is  crucial. 
In  this  case,  the  arrangement  is  such  that  many  of  the  levels  are  pushed  above  tf.  A  similar 
interaction  in  molecular  systems  will  bring  about  a  gain  in  the  occupation  of  one  of  the 
levels  at  the  expense  of  the  other,  but  here,  because  Cf  is  fixed,  both  metal  xz  and  CO  2x* 
levels  lose.  An  accompanying  diagram  could  be  drawn  for  the  CO  5 a  interaction  with  the 
bottom  of  the  xs  band  (5b).  However,  most  of  the  levels  will  lie  much  below  6f,  so  that 
even  an  interaction  of  similar  magnitude  will  not  alter  the  5 a  occupation  to  the  extent  of 
the  2x\ 

Should  this  interaction,  by  which  S  modifies  the  ability  of  some  surface-localized  metal 
states  to  bond  to  the  2x*  of  CO,  be  called  a  direct  one?  We  would  prefer  to  refer  to  it  as 
“through-bond  coupling” ,  a  concept  that  has  served  well  in  organic99*'*  and  inorganic94*  * 
chemistry.  Seemingly  localised  orbitals  (lone  pairs  in  two  parts  of  an  organic  molecule,  or 
two  metal  d  orbitals  in  a  binuclear  complex  connected  by  one  or  many-atom  bridges)  in 
two  or  more  parts  of  a  molecule  can  mix,  interact,  and  split  in  energy  as  a  result  of  mixing 
with  orbitals  of  intervening  atoms  or  groupings  of  atoms.  This  is  what  happens  here. 

That  the  5 a  interaction  is  modified  at  all  can  best  be  followed  indirectly  through  the 


Crystal  Orbital  Overlap  Population  (COOP)  curve3*  of  the  bridged  nickel  atoms.  The 
COOP  is  created  by  weighing  the  DOS  by  its  contribution  to  the  Ni-Ni  o.p.  In  Figure  5, 
these  are  compared  for  the  clean  and  sulfided  systems.  From  5*  and  5b,  we  know  that 
some  Ni-Ni  bonding  states  (near  the  bottom  of  the  d  band)  can  interact  with  the  5<r,  and 
that  some  antibonding  ones  interact  with  the  2x*.  Indeed,  the  resonances  with  the  CO 
are  picked  up  and  the  shift  from  bonding  to  antibonding  through  the  d  block  occurs  as 
expected.  However,  a  small  bonding  region  (circled  in  Figure  5)  appears  above  Cf  of  the 
sulfided  slab.  Assuming  conservation  of  the  absolute  number  of  bonding  and  antibonding 
states,  this  packet  must  have  been  pushed  up  from  the  d  block  by  the  sulfur  and  again  by 
the  CO  5<r.  Surface  levels  that  should  have  been  filled  by  o-donation  remain  empty.  The  a 
contribution  to  the  chemisorptive  bond  must  be  reduced,  regardless  of  the  stability  of  the 
5 a  occupation.  Since  the  number  of  affected  levels  is  not  substantial,  the  2 x*  occupation 
remains  the  overriding  factor. 

Figure  5  here 

Can  analogous  arguments  be  used  to  explain  the  anti-Blyholder  behavior  of  the  4-fold 
tetra  system?  The  charge  densities  computed  for  the  coadsorbed  slab  lie  much  closer  to 
the  single  species  values  than  do  their  bridging  counterparts.  A  few  are  listed  in  Table  6. 
That  the  magnitudes  are  small  should  not  concern  us  unduly  as  the  trends  tend  to  be 
more  reliable  in  similar  extended  Huckel  calculations.  Moreover,  the  trends  are  consistent 
among  themselves.  For  example  the  2 x*  occupation  remains  directly  proportional  to  the 


Ni-Ni.  with  S  &  CO 
Ni-Ni,  with  CO 


The  4-fold  tetra  system  differs  from  the  others  in  several  respects.  It  is  unique  among 
the  coadsorbed  systems  in  that  S  and  CO  are  each  able  to  withdraw  more  electron  density 
from  the  substrate  than  if  absorbed  individually.  It  is  also  the  only  system  for  which  the 
electronic  shifts  in  the  surface  levels  caused  by  the  coadsorbates  are  nearly  additive.  The 
sum  of  the  charge  transfer  between  the  bare  and  single  adsorbate  Ni(100)  surfaces  (columns 
2  plus  4  of  Table  6)  is  nearly  equivalent  to  the  shifts  between  the  bare  and  coadsorbed  values 
(column  3).  The  largest  difference  is  in  the  degenerate  xx  and  ys  orbitals.  Individually, 
S  and  CO  remove  a  total  of  0.122e~  from  each,  but  if  coadsorbed,  0.014e~  more.  As  a 
comparison,  the  greatest  discrepancy  for  the  bridged  CO  systems  is  nearly  five  times  the 
value  (in  xJ).  In  addition,  the  coadsorbates  work  on  entirely  different  parts  of  the  surface. 
This  can  be  seen  by  comparing  the  relative  charge  transfers  between  coadsorbed  and  single 
species  slabs  to  coadsorbed  and  bare  slabs  (columns  2-5  in  Table  6).  The  effect  of  CO  on 
the  clean  or  sulfided  substrate  is  greatest  at  the  surface  xy,  and  smallest  at  the  xz  and 
yz.  The  reverse  is  observed  for  the  sulfur  modification  of  the  clean  or  CO  preadsorbed 
slabs.  Thus  the  adsorption  characteristics  of  either  species  are  nearly  independent  of  the 
pertubation  of  surface  levels  caused  by  the  other.  The  same  information  can  be  obtained 
by  overlaying  the  appropriate  projected  DOS  curves  and  integrals. 


Table  6.  Electron  Density  Changes^,  4-Fold  Tetra  CO 


Ni(100)-p(2x2)CO 

Ni(100)-p(2x2)S+p(2x2)CO 

Atom(s)  or 

relative  to 

relative  to 

Orbitals 

Ni(100)  or  CO 

Ni(100)-p(2x2)CO 

Ni(100) 

Ni(100)-p(2x2)S 

CO  total 

+0.777 

+0.064 

- 

— 

CO  5 a 

-0.368 

+0.009 

— 

— 

CO  2jt* 

+0.694 

+0.023 

— 

— 

S  total 

— 

— 

— 

+0.015 

S  x,y 

— 

— 

— 

Surface  atom 

-0.283 

-0.338 

-0.621 

-0.266 

zl 

-0.043 

-0.071 

-0.114 

-0.044 

xy 

-0.133 

-0.023 

-0.154 

-0.132 

xz,yz 

+0.003 

-0.129 

-0.126 

-0.011 

‘Net  electron  gain  if  positive,  loss  if  negative. 


The  answer  may  lie  in  the  middle  layer,  rather  than  the  surface.  Adsorption  of  either 
S,  CO  or  both  will  increase  the  electron  density  of  the  bulk  layer  relative  to  that  of  the 
bare  slab.  The  effect  is  approximately  additive  (i.e.,  the  sum  of  the  single  species  = 
coadsorbed  system)  in  every  case  save  the  4-fold  tetra  site.  Here,  the  electron  transfer  is 
slightly  more  than  half  as  much  as  expected.  A  comparison  of  the  bulk  atom  projected 
DOS  curves  uncovers  a  striking  variation  in  the  behavior  of  the  xs  and  ys  on  the  bulk 
atom  lying  immediately  below  4-fold  coordinated  CO.  These  levels  do  not  interact  with 
2x*  of  either  the  single  species  or  the  4-fold  ortho  coadsorbed  system.  However,  nearly 
10%  of  the  xs  and  ys  states  are  pushed  up  into  the  4-fold  tetra  2 x*  peak.  In  simple  terms, 
an  additional  metal  to  CO  back  bonding  interaction  has  been  “turned  on”,  thus  the  2  ir* 
occupation  increases.  Why  is  this  particular  interaction  available  at  this  site?  It  is  not 
due  simply  to  the  tetragonal  symmetry,  since  the  bulk  xs  and  ys  do  not  interact  with  2** 
in  a  Ni(100)-c(2x2)CO  system.  The  later  has  the  same  total  adsorbate  coverage  as  the 
coadsorbed  systems;  each  S  is  simply  replaced  by  an  additional  CO.  We  can  speculate  that 
the  effect  of  the  S  on  the  bulk  atom  below  the  4-fotd  tetra  CO  is  stong  enough  to  modify 
the  levels  in  some  way,  but  too  weak  to  constrain  their  interactions. 

Clearly,  the  mechanisms  for  interadsorbate  interactions  are  severely  limited  in  this  sys¬ 
tem.  Not  only  does  the  distance  preclude  a  direct  effect,  but  adsorbate-surface  interactions 
are  segregated  such  that  the  potential  strength  of  an  indirect,  substrate  mediated  effect  is 
minimized  as  well.  At  this  point,  a  good  explanation  of  the  trends  is  still  obscure,  and  we 
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hesitate  even  to  speculate  on  its  origins. 

The  effect  of  sulfur  on  CO  adsorption  appears  to  be  determined  by  the  interadsorbate 
distance.  At  short  separations  the  S  and  CO  interaction  is  repulsive  —  we  see  site  blockage. 
At  intermediate  separations,  for  instance  in  the  energetically  favorable  geometry  where  S  is 
in  a  four-fold  hollow  and  CO  is  bridging,  we  see  the  workings  of  the  usual  Blyholder  model. 
The  2x*  of  CO  is  populated  less,  with  resultant  weakening  of  C-Ni  and  strengthening 
of  C-0  bonds.  We  trace  this  effect  to  the  modification  by  S  of  the  ability  of  specific 
surface  orbitals  (xs  and  ys)  to  backbond  to  2 *•*,  or  through-bond  coupling.  At  longer 
separations,  we  see  a  small  reverse  effect  (C-0  weaker,  C-Ni  stronger)  which  has  some 
experimental  support,  but  it  is  not  easy  to  rationalise.  In  general,  we  think  that  the 
simple  electronegativity  and  electron  transfer  explanation  of  coadsorption  effects  needs  to 
be  supplemented  by  a  detailed  orbital  analysis  in  terms  of  through-space  and  through-bond 
interactions. 
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Appendix 

The  calculations  were  performed  using  the  tight  binding  extended  Huckel  method2  u,k. 
The  Hu’s  for  Ni  had  been  previously  determined2*  by  charge  iteration  on  the  bulk  metal 
using  Gray’s  equations37.  The  Hu’s  for  C  and  O  had  also  been  determined  previously20  by 
three  cycle  iteration  on  CO  adsorbed  onto  a  four  layer  Fe(110)  slab.  A,  B,  and  C  iteration 
parameters  are  from  reference  38.  All  parameters,  including  the  uniterated  S  parameters, 
are  listed  in  Table  7. 

Table  7  here 

For  both  the  Ni(l00)-p(2x2)S  and  c(2x2)S  systems,  the  Ni-S  bond  length  is  reported 
to  be  2.19A30.  The  Ni(100)-p(2x2)CO  geometry  is  d(Ni-C)=1.80A  and  d(C-0)=1.15Al2c. 
Adsorbates  were  placed  on  one  side  of  the  three  layer  slab.  Ten  and  16  k  point  sets  were 
used  respectively  for  systems  of  tetragonal  and  orthorhombic  symmetry.  Convergence  was 
checked  for  10,  15  and  28  k  point  sets  on  the  Ni(100)-p(2x2)S  system.  The  maximum 
variations  among  these  sets  were  0.007eV  for  the  Fermi  energy,  O.OOleV  for  the  total 
energy,  0.010  for  orbital  electron  densities  and  0.001  for  overlap  populations. 


Orbital 

Hu(eV) 

ft 

ft 

<Y 

<Y 

Ni  4s 

— 7.S 

2.1 

4p 

-3.7 

2.1 

3d 

-9.9 

5.75 

2.00 

0.5683 

0.6292 

C  2s 

-18.2 

1.63 

2p 

-9.5 

1.63 

0  2s 

-29.6 

2.27 

2p 

-13.6 

2.27 

S  3s 

-20.0 

1.817 

3p 


-13.3 


1.817 
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Figure  1.  In  the  left  panel,  the  bare  Ni(lOO)  DOS.  The  valence  orbitals  of  CO  and  S  are 
drawn  in  the  middle  and  right  panels  respectively. 

Figure  2.  The  solid  line  indicates  the  Ni(100)-p(2x2)S  total  DOS.  The  dotted  and  dashed 
lines  are  the  integrals  of  the  surface  xx  +  yz  and  S  x  +  y  projected  DOS  respectively. 

The  S  unsupported  square  net  DOS  is  represented  by  the  median  energy  bars  on  the 
right. 

Figure  3.  The  bridging  Ni(100)-p(2x2)S+ p(2x2)CO  total  DOS;  major  peaks  are  labeled. 

The  bars  on  the  right  indicated  the  median  energies  of  the  CO  and  S  unsupported 
nets. 

Figure  4.  the  dashed,  dot-dashed  and  dotted  lines  are  respectively  the  integrals  of  the  xz 
projected  DOS  of  Ni(100),  Ni(lOO)-p(2x2)CO  and  Ni(100)-p(2x2)S+p(2x2)CO.  The 
solid  curve  is  the  total  DOS  of  the  latter.  The  median  energies  of  the  5 a  and  2rr* 
coadsorbed  levels  are  indicated  on  the  right. 

Figure  5.  The  bridged  Ni-Ni  COOP’s  for  the  Ni(lOO)-p(2x2)CO  and  Ni(100)-p(2x2)S+p(2x2)CO 
are  indicated  by  solid  and  dotted  lines  respectively. 


